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Abstract 

A Rayleigh scattering diagnostic for high speed 
flows is described for the simultaneous, instantaneous 
measurement of gas temperature and velocity at a 
number (up to about one hundred) of locations in a 
plane illuminated by an injection-seeded, frequency 
doubled Nd:YAG laser. Molecular Rayleigh scattered 
light is collected and passed through a planar mirror 
Fabry-Perot interferometer. The resulting image is 
analyzed to determine the gas temperature and bulk 
velocity at each of the regions. The Cramer Rao lower 
bound for measurement uncertainty is calculated. 
Experimental data is presented for a free jet and for 
preliminary measurements in the Lewis 4 inch by 10 
inch supersonic wind tunnel. 

Introduction 

Quantitative techniques for measurement of fluid 
flow properties in aerospace test facilities are needed to 
obtain performance data, to validate new computer 
codes, and to establish facility flow conditions. In many 
experiments, conventional probe measurements are not 
adequate because they perturb the flow under study. 
This is a particularly severe problem in high speed 
flows in small wind tunnels and in other confined flows 
such as engine components. 

Techniques to measure velocity based on Mie 
scattering, such as LDV and PIV, require micron size 
particles to be entrained in the flow. These particles 
may not always accurately track the flow under study, 
and thus the velocities obtained may not represent the 
actual fluid velocity. Particle lag errors can be large 
when the fluid undergoes high accelerations, as in flows 
with shocks. Also, in some applications it is either not 
feasible or not desirable to introduce foreign material 
into the flow. For example, it is difficult to seed large 
wind tunnels without contaminating the tunnel. Very 
high temperature flows, such as in rocket combustion 
chambers, also make seeding difficult. For these 
reasons, techniques that are based on molecular 
scattering are of interest. 
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A variety of molecular scattering diagnostic 
techniques have been developed for the measurement of 
gas density, temperature, and velocity. Included are 
Rayleigh and spontaneous Raman scattering 12 , laser 
induced fluorescence 3 , laser flow tagging 4 , and transient 
grating spectroscopy 5 . 

Rayleigh scattering based diagnostics, the subject 
of this paper, offer several important features when 
compared to other molecular scattering methods. 
Rayleigh scattering is a nonresonant process, so any 
wavelength laser may be used. This generally means a 
lower cost, simpler laser system. And, Rayleigh 
scattering is much stronger (typically by two orders of 
magnitude) than spontaneous Raman scattering. On the 
other hand, the resonant processes (such as laser 
induced fluorescence) generally are much stronger, but 
they require more complex, more expensive tunable 
laser systems. Because Rayleigh scattering is not 
species specific, the gas composition must be known for 
density and temperature measurements. This, however, 
is not a factor for an air flow diagnostic, as is being 
addressed in this paper. 

The simplest application of Rayleigh scattering is 
to measure the total amount of scattered light, which is 
proportional to the gas density. When the pressure is 
known, this technique can be used to obtain temperature 
data. This approach is susceptible to interference from 
stray laser light and to Mie scattering from particles 6 . 
(However, use of two wavelengths can be used to 
reduce the interference from stray laser light 7 .) More 
information is available from the spectrum of the 
Rayleigh scattered light. The spectrum is directly 
related to the velocity distribution of the molecules and 
thus contains information about the thermodynamic 
parameters, including temperature and bulk velocity. 

The spectral width of Rayleigh scattered light is a 
result of Doppler broadening caused by the motion of 
the gas molecules. Since typical molecular speeds are 
on the order of the speed of sound, the line width of the 
Rayleigh line is only about 1 GHz for air and 90° 
scattering. A high resolution instrument is thus required 
to measure the Rayleigh spectrum. Narrow band 
absorption filters and Michelson and Fabry-Perot 
interferometers have been used for this work. 

Measurement of the Rayleigh scattering spectrum 
with a scanning Fabry-Perot interferometer has been 
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applied to a variety of flow and combustion 
experiments, such as determining time-average 
temperature and/or velocity at a point in a hydrogen air 
flame 8,9 , in a hydrogen oxygen rocket exhaust 
plume 1011 , and in a wind tunnel test of an ASTVOL 
model 12 . 

The Fabry- Perot interferometer has also been used 
in a static imaging mode (rather than a scanning mode) 
to measure velocity and temperature in the atmosphere 
by observing the emission spectrum of atomic oxygen 13 . 
This imaging method has been used with Rayleigh 
scattering for one dimensional and planar 
temperature 14 15 and velocity measurements 16 . 

Molecular absorption filters are used for frequency 
discrimination in the Filtered Rayleigh Scattering 
technique 17 , which was developed to measure density, 
temperature, and velocity over a planar region. An 
angularly resolved filtered Rayleigh scattering 
technique 18 is also being developed for air mass flux 
measurements and other high speed flow studies. 

In this paper, our previously described technique 16 
based on a CW laser and Fabry-Perot interferometer to 
determine time-average velocity in a supersonic free jet 
over a 2D region is extended to instantaneous 
measurement of velocity and temperature. The beam 
from an injection- seeded, frequency doubled Nd:YAG 
laser is formed into a sheet that illuminates a planar 
region in the flow. Molecular Rayleigh scattering is 
collected from the flow and imaged through a planar 
mirror Fabry-Perot interferometer onto a gated, 
intensified CCD detector. Analysis of the recorded 
interference pattern provides estimates of the gas 
temperature and one component of the velocity at about 
one hundred localized regions in the imaged area. The 
use of a pulsed laser for this work allows the 
determination of the instantaneous velocity field which 
can be used to calculate turbulence parameters and 



spatial correlations. In addition, the pulsed laser 
permits easy identification of particle images, which can 
be eliminated in the data processing. An analysis is 
performed to calculate the Cramer Rao lower bounds 
for the measurement uncertainties as a function of the 
flow parameters and the measurement system 
parameters. Experimental results are presented for a 
free jet at high subsonic velocities, and some 
preliminary results are shown for the Lewis 4 inch by 10 
inch supersonic wind tunnel. 


Theory 

Ravleigh scattering 

The most basic application of Rayleigh scattering 
as a diagnostic is to determine gas density by measuring 
the total amount of Rayleigh scattered light (fig. 1). 
This does not require resolution of the spectrum. 

However, because the Rayleigh scattering cross 
section is species dependent, the gas composition must 
be known. If the gas pressure is known, the density can 
be used to determine temperature. The number of 
Rayleigh scattering photons N R collected by an optical 
system with solid collection angle Q is 




E.nL x XQ( d<j\ . 2 
—2 — * — — sim x 

he U n) 
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where E fl is the laser energy, A is the wavelength, n the 
gas number density, L x the length along the beam of the 
scattering volume, h Planck's constant, c the velocity of 
light, da/dH the differential scattering cross section, 
and x the angle between the electric field vector of the 
(linearly polarized) incident light and the direction of 
the scattered light. For example, about 500 photons are 
collected by an f/8 optical system from a 1 mm 3 
volume of nitrogen at NTP irradiated by a 1 mJ, 532 nm 
laser pulse. 

Resolving the spectrum (fig. 2) of the Rayleigh 
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Fig. I - Light scattering experiment 


Fig. 2 - Rayleigh scattering spectrum 
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scattering gives information on gas velocity and 
temperature. The gas velocity can be determined from 
the shift in the peak of the spectrum relative to the 
frequency of the incident laser beam. This shift is 
proportional to one component of the mean velocity, in 
a manner analogous to the Doppler shift encountered in 
laser anemometry. Additional velocity components can 
be obtained by changing either the direction of the 
incident beam or the observation direction of the 
scattered light. Since velocity is determined from only 
the shift of the spectral peak, this measurement is not 
affected by the gas composition. Also, scattering from 
particles entrained in the flow will have a spectral peak 
at the same location as the Rayleigh scattering peak. 

Gas temperature is related to the width of the 
spectral peak. The shape of the Rayleigh spectrum is a 
function of the scattering angle and gas density. The 
simplest example is a single component, low density gas 
with a Maxwellian velocity distribution. In this case, the 
scattering occurs independently at each molecule, and 
the Rayleigh scattering spectrum is a simple Gaussian 


S(f)df = 



'27T(/-/„)-K-u' 
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where, as shown in figure l, K = ks-ko is the interaction 
wave vector with magnitude (wave number) K - 
(4/A)sin(0 5 /2), k* and ks are the wave vectors of the 
incident and scattered light, 0 S is the scattering angle, a 
= (2KT/m) m is the most probable molecular speed (with 
ic being Boltzmann's constant, m the molecular mass, T 
the gas temperature), and u is the mean gas velocity. 

For higher density gases (where the molecular 
mean free path becomes comparable to or smaller than 
the scattering interaction wavelength 2nlK), the 
spectrum is no longer Gaussian in shape. The ratio of 
the interaction wavelength to the mean free path is 



Fig. 3 - Rayleigh scattering spectra for three values of y parameter 
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Fig. 4 - Fabry-Perot interferometer 


normally expressed as the parameter y = p!r\Ka, where 
p is the gas pressure and T] is the shear viscosity. The 
Gaussian regime (low density) corresponds to y « 1- 
For v » 1 (high density), the scattering spectrum is 
strongly influenced by collective effects and is 
characterized by a central peak and two sidebands. The 
sidebands can be thought of as being caused by 
scattering from thermally excited acoustic waves. A 
continuum theory 19 can be used to model the spectrum 
here. The spectrum in the transition regime, where y ~ 
I, requires a more detailed kinetic theory. We use the 
Tenti S6 model 20 for this work. Typical spectra for the 
three regimes are shown in Figure 3. 


Fabrv-Perot interferometer 

The Fabry-Perot interferometer (as shown in fig. 
4) consists of two partially transmitting planar mirrors. 
Multiple reflections between the mirrors result in a 
transmission function (defined as the fraction of light 
transmitted by Fabry-Perot for a monochromatic source) 
given by 21 




1 + Fsin 


)( 2k f {id a cosfl r 


(3) 


Here, f=cfK> is the optical frequency, ju is the refractive 
index of the medium in the Fabry-Perot cavity, d 0 is the 
Fabry-Perot mirror spacing, 0 r is the angle between the 
ray and the optic axis, and F = 1 /(sin“(Tc/2 A^ e) where Ne 
is the effective finesse. The image of a monochromatic 
extended source located in the object plane consists of a 
series of concentric rings, such as shown in figure 8. 

Ravleigh spectral measurement using Fabry-Perot 
interferometer 

Consider a laser sheet of energy E 0 that uniformly 
illuminates a region in the object plane (fig. 5). The 
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Fig. 5 - 2D scattering experiment using imaging mode Fabry-Perot 
interferometer 


a molecular absorption filter. Finally, S R (f&) is the 
spectrum of the Rayleigh scattered light, which is a 
function of the flow parameters and the optical 
parameters. The integration over the solid angle is 
necessary to account for the range of K vectors over 
the light collection aperture. For small solid collection 
angles or low velocities, this effect can be neglected 10 
and the integration in equation 4 over £2 can be 
eliminated. For the work presented here, the effect of 
aperture broadening was calculated to be on the order 
1 % or less and was neglected. Spectral broadening can 
also be caused by sufficiently strong turbulence, but it 
was not considered in this work. Velocity measurements 
are not sensitive to aperture broadening since only the 
frequency of the spectral peak is used, and the spectrum 
is symmetrically broadened. 


height of the sheet is assumed to be equal to the height 
of the area viewed by the array detector. The array has 

N r rows (measured perpendicular to the beam 
propagation direction) of square pixels of size L x by L x . 
The expected number of detected photons for the q th 
pixel can be written 

... 

+ < A ' Wg >S(f-f n )]/(/) d/d4rfQ+ B q 


where 


< N Rq >dQ = — — ■ — 


N r he M 


— |</n 
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is the number of photons scattered (Rayleigh scattering) 
into solid angle dQ that would be detected by the q * 
pixel if the Fabry-Perot interferometer were not present. 
Similarly, <N w f>dQ is the number of photons due to 
laser light scattered from windows and walls and is at 
the laser frequency / 0 =c/A* f . Broadband background 
light, detector dark current, and readout noise are 
represented by B q . Although equation 4 does not 
include scattering from particles, it could be easily 
incorporated. In equation 5, M is magnification (i.e., 
size of image of the CCD in the measurement region is 
(N r xN c )L x /M where N c is the number of columns), and e 
represents the overall detection efficiency (including 
detector quantum efficiency and system losses). The 
instrument function 1(f) is made up of the Fabry-Perot 
transmission function /fp(/) and any other filters, such as 


Data processing 

The Rayleigh scattering images are processed 
using a nonlinear least-squares fit to a model function 
represented by equation 4. The y parameter for this 
flow condition and optical configuration is on the order 
of one, so the Tenti S6 model is used for calculation of 
the Rayleigh scattering spectrum. The necessary 
thermodynamic parameters are calculated using the 
FLUID program 22 . 


The data reduction procedure is as follows: a 
number of 30x30 pixel subregions in the recorded 
image are defined on the fringes (a typical set of 
subregions is shown in figure 6). A reference image (of 
elastically scattered laser light) is first analyzed to 
determine the interferometer finesse and phase of the 
unshifted laser light at each subregion. (The measured 
finesse includes the effect of the finite line width of the 
Nd:YAG laser.) The data obtained at the set of 
subregions are fit to a 2D polynomial. The image of the 



Fig. 6 - Set of subregions used for data processing (actual subregion 
size is 30x30 pixels, although shown here as 15x15 pixels for 
clarity). 
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Rayleigh scattered light is then analyzed at each 
subregion using the finesse and phase obtained from the 
reference image. The velocity is calculated from the 
change of phase and the temperature from the spectral 
broadening. Another factor that must be accounted for 
is the broadening caused by the multichannel plate 
intensifier; the full width half maximum spread of a 
point source is about 3 pixels. This effect is included 
by convolving a 2D Gaussian function with the image. 
As a result of the extensive calculations required for the 
full evaluation of the model function for the iterative 
fitting procedure, the time is quite long for the 
processing an entire image — about 2 minutes on a 90 
MHz Pentium PC for each 30x30 pixel subregion. The 
result of this procedure is a value for the velocity 
component and the gas temperature at each subregion. 
For presentation purposes, the sets of velocities and 
temperatures can be converted to contour plots. It 
should be reiterated that velocity and temperature data 
are only obtained for subregions located on the fringes. 
Data are obtained at a large number of locations, but no 
data are obtained in the regions between fringes (see 
fig. 8). To achieve more complete coverage, additional 
images can be obtained at different Fabry-Perot mirror 
spacings, which give fringes at different radii. (The 
data are then not all obtained simultaneously.) In any 
case, the results obtained with this procedure for each 
subregion represent average values of velocity and 
temperature over the subregion. (The size of the 
subregions for this work are about 0.4 mm square.) 



( 6 ) 


where N R is the total number of detected photons due to 
Rayleigh scattering (eq. 1 or eq. 5). Stray laser light 
and background noise are neglected. 

As an example, consider nitrogen at NTP (T=293 
K, p = 1 atm., n = 2.5x1 0 25 m' 3 ), a single 1 J pulse from 
a 532 nm laser, do/dQ = 6. lx 10 32 m 2 /sr, an array 
detector having 384 rows of 23 pm square pixels, a 1 : 1 
imaging system having f/8 collection optics (Q. = 
0.0123 sr), and a 1% overall detection efficiency. 
Using equation 5, the expected number of detected 
photons from a single pixel area is <Nr q > = 30. Thus, 
for a 30x30 pixel region, the number of expected counts 
would be 27,000, and the lower bounds on the precision 
of density, temperature, and velocity measurements are 
0.6 %, 2.5 K, and 1.8 m/sec, respectively. It must be 
emphasized that these uncertainties are based on the 
assumption of ideal spectrum analysis. Actual 
measurements based on some type of interferometer or 
absorption medium will, of course, have a higher lower 
bound for the uncertainties. It was previously shown 
that for the case of no wall scatter, the lower bound for 
the Fabry-Perot based system is about a factor of 5 


Uncertainty analysis 

For light scattering experiments, the uncertainty 
caused by photon shot noise can readily be evaluated. 
If this uncertainty is the dominant source of 
measurement error, as it is in many experiments, then 
this analysis can be used to estimate the measurement 
error, and it can be used as a tool for designing the 
experiment to minimize the error. (On the other hand, if 
the shot noise induced uncertainty is small compared to 
other error sources, this analysis comprises only one 
part of the overall uncertainty analysis.) In any case, 
the shot noise sets the lower limit for measurement 
uncertainties. The theoretical least possible uncertainty 
in the measurements is given by the Cramer Rao lower 
bound 23 . For an ideal spectrum analyzer, the lower 
bounds for the uncertainty in the density, measured 
velocity component, and temperature for a Gaussian 
Rayleigh spectrum are 16 



Fig 7 - Optical setup for measuring axial velocity component in free 
jet 
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larger than for an ideal spectrum analyzer for velocity 
and about a factor of 10 larger for temperature. 

Experiments 

Free jet 

Initial laboratory measurements were made in a 
small free jet to validate the technique (fig. 7). The exit 
diameter was 9.5 mm. The air supply to the nozzle was 
filtered to remove particulates with diameters larger 
than 0.2 pm. The beam from an injection seeded, 
frequency doubled, Nd:YAG laser with 0.9 J pulse 
energy was formed into a sheet that crossed the flow at 
45° as shown in the figure. The sheet thickness was 
about 0.2 mm and the height was about 7 mm. A 
commercially available Fabry-Perot interferometer was 
used. The mirrors (70 mm diameter, A/200 flatness, and 
90% reflectivity) were set at a 14.96 mm spacing for the 
measurements reported here, which corresponded to 
free spectral range (FSR) of about 10 GHz. Rayleigh 
scattered light was collected and collimated with lens 
LI (250 mm focal length lens). The scattering angle 
was 90°. After passing through the pellicle beam 
splitter and Fabry-Perot interferometer, the light was 



Fig. 8 - Reference image for first free jet image 



Fig. 9 - First Rayleigh scattering image for free jet; note particle 
images in upper right 



Fig. 10 - Velocity contour plot for first free jet image; contours in 
m/sec 

focused onto the detector (multichannel plate 
intensified, cooled CCD array with 576x384 23 (im 
square pixels) by lens L2 (300 mm focal length 35 mm 
camera lens with 1.4X teleconverter). The images 
digitized with a 14 bit A/D converter and transferred to 
a laboratory computer for storage and analysis. 

A beamsplitter was used to direct a small amount 
of light from the laser onto an opal glass diffuser to 
serve as a reference image. A condenser lens 
collimated the light from the diffuser, which was then 
directed into the Fabry-Perot by an uncoated pellicle. 
Remotely controlled shutters were used to block either 
the main beam (for recording the reference image) or 
the reference beam (for recording the Rayleigh image). 

The measurement plane was intersected the by jet 
axis 12 mm downstream from the nozzle exit plane. 
Data for a single flow condition are presented here; the 
axial velocity and static temperature calculated from the 
isentropic flow equations were 255 m/sec and 280 K. 
Three representative images are shown. Each was taken 
with a single laser pulse. The Fabry-Perot mirror 
spacing was different for the three images, so the 
fringes are at different radii. The first Rayleigh image 
along with its associated reference image are shown in 
figures 8 and 9. This image was analyzed at the 90 
subregions shown in fig. 6. The flow direction was 
assumed to be along the jet axis for the purpose of 
converting the measured velocity component (which is 
along K) to the velocity magnitude. The y parameter 
was about 0.9, so the Tenti S6 model for the Rayleigh 
spectrum was used. Velocity data were converted to the 
contour plot shown figure 10. The velocity 
measurements near the center of the image were close to 
the isentropic flow values. Temperature data are not 
shown for this image because of inadequate detector 
resolution at the outer fringes. 
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Fig. 1 1 - Second free jet Rayleigh scattering image 


Fig. 13 - Third free jet Rayleigh scattering image 



Fig. 12 - Subregions analyzed for second free jet image (fig. 1 1 ) 

The second and third Rayleigh images were 
analyzed at a smaller number of subregions located on 
the inner fringes. The second image was analyzed at 30 
subregions (fig. 1 1). The mean velocity was 250 m/sec, 
with standard deviation of 1 1 m/sec. The mean 
temperature was 270 K, with standard deviation of 23 
K. 

The third image was analyzed at 19 subregions 
(fig. 12). The mean velocity was 252 m/sec, with 
standard deviation of 1 1 m/sec. The mean temperature 
was 267 K, with standard deviation of 16 K. 

For this optical geometry and flow, the Cramer 
Rao lower bounds (for an ideal spectrum analyzer) for 
velocity and temperature were about 0.4 m/sec and 0.6 
K. The measurement precision based on the observed 
variance in the data for the second and third images was 
about 10 m/sec and 20 K. The variance was greater for 
the first image with more subregions covering a larger 
part of the image field. The might be caused by the 
limited spatial resolution of the intensified CCD. Good 
agreement between the measured mean velocity and the 
velocity calculated from the isentropic flow equations. 


Fig. 14 - Subregion analyzed for third free jet image (fig. 13) 

However, the mean measured temperature was 
significantly less than the isentropic temperature, 
particularly in the first image with 90 subregions where 
it was 50 K low. The the other two images, with fewer 
subregions, were low by about 10 K. Note that no 
adjustable parameters were used in obtaining these 
velocities and temperatures. Measurements of a known 
temperature gas could be used to obtain calibration 
parameters that could be used to increase the accuracy 
of the temperature measurements. 

Supersonic wind tunnel 

The Rayleigh scattering system was recently 
installed in the Lewis 4 inch by 10 inch supersonic wind 
tunnel, and some preliminary data were obtained. The 
optical system was identical to the layout described in 
the above section for the free jet measurements, except 
for the optics used to introduce the beam into the tunnel. 
The tunnel was equipped with a Mach 2 nozzle, and the 
air supply was passed through a 0.2 pm filter to remove 
particles. For the first tests, the laser beam was 
transmitted through a window into the plenum chamber 
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Fig. I? - Rayleigh scattering image for supersonic wind tunnel 



Fig. 16 - Subregions (27) analyzed for supersonic wind tunnel image 



Fig. 17 - Temperature contours for supersonic wind tunnel (fig. 15) 

and propagated downstream through the nozzle and the 
test region. Light scattered at 90° was collected through 
a window in the tunnel sidewall. A typical fringe 
pattern (for a single laser pulse) is shown in figure 15. 


For this measurement, the laser beam was focused to 
about 2 mm diameter. As can be seen in the image, the 
amount of stray laser light detected was small, and the 
Filter effectively removed particles in the air supply. 

One problem in this test environment was the high 
level of acoustic noise, which caused vibration of the 
optics. This resulted in shot-to-shot variation in the 
phase of the reference images. As a result of this 
variation, the velocity could not be obtained from the 
phase shift between subsequent reference and Rayleigh 
images. Instead, a value for the reference phase was 
obtained from a part of the Rayleigh image dominated 
by stray, unshifted laser light (e.g., near the bottom of 
the image in fig. 13). Temperature data, since they are 
obtained from the broadening rather than the location of 
the spectral peaks, are not as sensitive to the shot-to- 
shot phase variation. A temperature contour plot based 
on 27 subregions (fig. 16) is shown in figure 17. The 
mean measured temperature was 202 K with a standard 
deveiation of 23 K. The estimated velocity was 390 
m/sec, which correspond to about Mach 1 .4 flow based 
on isentropic flow from the plenum. Acoustic shielding 
material will be used in future test to isolate the laser 
and optics. Another approach, which was used for 
these preliminary measurements, would be to obtain the 
reference frequency information from the same image 
used for the Rayleigh light. 

Concluding Remarks 

Molecular Rayleigh scattering measurements made 
in a free jet and in a small supersonic wind tunnel 
demonstrated the feasibility of using an injection-seeded 
Nd:YAG laser and Fabry-Perot interferometer to obtain 
a large number (on the order of one hundred) 
simultaneous, instantaneous measurements of gas 
temperature and velocity in a planar region of flow 
fields. Good agreement was obtained in free jet 
measurements for both temperature and velocity. 
Several issues, however, must be addressed in the 
application of this technique in the harsh environment 
of test facilities. 

Stray laser scattering must be carefully controlled. 
In many experiments the stray laser light can be 
managed with apertures and light traps. In high speed 
flows (supersonic), where the spectrum is shifted from 
the laser frequency, a molecular absorption filter can be 
used to reject the stray light. In addition, the laser and 
the Fabry-Perot interferometer must be isolated from 
high noise levels. This can be done with acoustic 
shielding and vibration isolation equipment, or the laser 
and interferometer can be remotely located and the light 
transmitted through relay optics or optical fibers. 
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